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Abstract. We show that toroidal magnetic field annihilation in the equatorial current sheet of the pulsar magne- 
tosphere is related to how fast poloidal magnetic field lines close as we move away from the light cylinder. This 
determines the radial reconnection electric field which directly accelerates particles in the neutral layer inside 
the equatorial current sheet. The efficiency of the poloidal magnetic field closure near the light cylinder may be 
measurable through the pulsar braking index. We argue that the lower the efficiency of pair formation, the higher 
the braking index. We also argue that the synchrotron radiation reaction in the neutral layer does not inhibit the 
accelerated particles from reaching the maximum energy of about ~ 10^^ eV that is available in the open pulsar 
magnetosphere. 
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1. Introduction 

After almost 40 years of pulsar research, the structure of 
the pulsar magnetosphere, the region extending from the 
neutron star surface to the wind termination shock, is still 
not clearly understood. The observed pulsed radiation 
(radio to gamma ray) is believed to originate in the mag- 
netosphere, and yet it is still not clear where and how it is 
produced (e.g. lArons 19831 Mushmov &: Harding 2004 ). 
The magnetosphere is also the source of a fast particle 
wind with Lorentz factor on the order of 10^ that 
terminates in a standing shock at a distance of about 
0.1 pc dBogovalov fc Aharonian 2000p . Several acceler- 
ation mechanisms have been proposed, ranging from 
ideal 



the toroidal magnetic field component is associated with 
an electric field that accelerates particles directly inside 
the current sheet. Although several authors have acknowl- 
edged the presence of this electric field (e.g. IKirk 20041 



MHD 

IVlahakis 2004D 
[Michel 1994( 
rect particle 
rent sheet 



models ( Contopoulos & Kazanas 2002 



pressure driving (|Coroniti 1990 



|Kirk fc Skjsraasen 20031 ), to di- 
acceleration in the equatorial cur- 
(Romanova, Chulsky & Lovelace 2005', 



Contopoulos 2007), but we still do not understand 



how this wind is produced. 

In recent years, there has been mounting evidence 
against ideal MHD wind acceleration, and there is no con- 
vincing self-consistent, ideal MHD wind solution available 
in the literature ( Bogovalov 1997[ ). Thus, we decided to 
direct our attention to models that contain regions where 
ideal MHD breaks down. The most promising of these re- 
gions is the equatorial current sheet, across which the main 
magnetic field component, the toroidal one, changes direc- 
tion and annihilates (Contopoulos 2007). Annihilation of 



Romanova, Chulsky fc Lovelace 2005| , we feel that a cer- 
tain important physical element has been neglected in 
the discussion, namely the global topology of the poloidal 
magnetic field. 

In |Contopoulos 2007[ we dealt with non-ideal MHD re- 
gions through the notion of the effective magnetic diffusiv- 
ity. In the present work, we give a simpler picture without 
mentioning the magnetic diffusivity in the equatorial cur- 
rent sheet. In § 2 we argue that toroidal magnetic field an- 
nihilation cannot be studied independently of the poloidal 
magnetic field. In § 3 we associate the pulsar braking in- 
dex to the rate at which poloidal magnetic-field lines close 
with distance in the vicinity of the light cylinder. Finally, 
in § 4 we discuss the issue of particle acceleration in the 
equatorial current sheet. 

2. The global topology of the pulsar 
magnetosphere 

Let us introduce a spherical system of coordinates (r, 0, (p) 
centered on the neutron star and aligned with its axis of 
rotation. The neutron star rotation delineates two distinct 
hemispheres, the upper and lower, separated by the rota- 
tion equator. The magnetic field also delineates two half 
spaces: the north, where the magnetic field points away 
from the star, and the south, where the magnetic field 
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points towards the star. In general, the rotation and mag- 
netic axes are at an angle ^ 0, and we may assume 
without any loss of generality that < 90°. In that case, 
Br is mostly positive in the upper ('north') hemisphere, 
whereas it is mostly negative in the lower ('south') hemi- 
sphere. 

In the case of an isolated pulsar, all magneto- 
spheric field lines naturally originate on the surface 
of the neutron star, since the source of the magneto- 
spheric field is the neutron star interior. Obviously, so- 
lutions also exist with field lines not connected to the 
star ( [Lovelace, Turner fc Romanova 2006] ) , and these may 
very well apply to neutron stars surrounded by disks 
with disk winds, but not to isolated neutron stars. The 
star holds the magnetic field from diffusing to infinity. 
The field close to the neutron star is dipolar, but as 
we move away from the stellar surface, the magnetic 
field topology is gradually distorted. The magnetic field 
is stretched outwards and is wound backwards because 
of the neutron star's rotation. The distortion depends 
on the magnetospheric plasma conductivity, on the in- 
clination angle between the rotation and magnetic axes, 
and, as we will see, on the neutron star spindown his- 
tory. The above characteristics can be found in sev- 
eral numerical solutions of simple limiting cases (e.g. 
Contopoulos, Kazanas fc Fendt 19991 [Spitkovsky 2006[ ). 



In the axisymmetric case (?? ~ 0), the open poloidal mag- 
netic flux in each magnetic hemisphere is given by 



Let us now discuss the significance of the magneto- 
spheric plasma conductivity, first considering the case of 
infinite plasma conductivity. Ideal MHD requires that 
plasma is confined to flowing along the magnetic field, 
so field lines may be visualized as corrotating with the 
neutron star at the stellar angular velocity fl. This does 
not mean that the magnetospheric plasma corrotates with 
the star. The plasma angular velocity is equal to O only 
wherever = (e.g. IBeskin 1997p . A direct corollary 
of this is that no field lines close outside the light cylin- 
der, which is the cylindrical distance tl = c/f2, where if a 
particle were to corrotate with the neutron star, it would 
do so at the speed of light. The argument is straightfor- 
ward: as we said, field lines from both stellar magnetic 
hemispheres are wound backwards by the neutron star's 
rotation; therefore, a point on every closed field line may 
be found where = 0. If a closed field line were to 
extend beyond the light cylinder, that point would lie 
outside the light cylinder, and according to our discus- 
sion above, the plasma velocity there would have to be 
at least equal to r^smO > c. The unavoidable conclu- 
sion is that the pulsar magnetosphere consists of closed 
field lines that corrotate with the neutron star inside the 
light cylinder and of field lines that cross the light cylinder 
and open up to infinity. It is important to keep in mind 
that nothing exceptional takes place at the light cylinder, 
because field lines cross it without kinks or discontinu- 
ities ( [Contopoulos, Kazanas fc Fendt 1999| . Several au- 
thors have shown that asymptotically (r>rL), the poloidal 
magnetic field structure becomes radial (monopole-like) . 



^'(r,6l) « - cosd) , 

and the magnetic field by 
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\Br{r,e)\^ Bl[ — ] , Be(r,0) «0 and 
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Bl sine 
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([Michel 19941 [Contopoulos, Kazanas fc Fendt 19991 ). 
Here, Bl = ^L/(27rr|^) is a typical value for the 
magnetic field at about the light cylinder distance. As 
|Bogovalov 1999[ shows in a very important and elegant 
paper, Eqs. [H are valid even in the case of an oblique 
(i.e. non-aligned) magnetic rotator. In that case, in the 
upper hemisphere, magnetic field lines mostly point 
away from the star {Br > 0, B^ < 0), whereas in the 
lower hemisphere they mostly point towards it {Br < 0, 
B^ > 0). The two regions are separated by a spiral 
current-sheet discontinuity described by Eq. 3 in Kirk & 
Lyubarski (2001), across which Br and B^f, change sign 
together. 

Let us next consider the more realistic case of high, 
but not infinite, plasma conductivity. It is natural to 
argue that ideal MHD is a realistic description of the 
largest part of the pulsar magnetosphere, except for re- 
stricted regions with high electric-current densities, where 



ideal MHD breaks down ( [Contopoulos 20071 ). I^i that 
case, field lines are not forbidden to close outside the 
light cylinder. In fact, several numerical solutions that 
try to simulate the ideal MHD magnetosphere show 
this different topology because of numerical diffusivity 
([Komissarov 20061 [McKinney 20061 [Spitkovsky 2006| . It 
is very important to note that field-line closure out- 
side the light cylinder may take place only if it is as- 
sociated with magnetic field energy dissipation (non- 
ideal MHD). In fact, several authors have proposed 
that this mechanism may be related to the accelera- 
tion of the pulsar wind, and this led to the so-called 
'striped' pulsar wind model ([Coroniti 19901 [Michel 19941 
[Lyubarsky fc Kirk 20011 [Kirk fc Skjaeraasen 20"03| . 

There is here, we believe, a subtle but very important 
point that has been overlooked by most previous studies 
of magnetic field dissipation in the current sheet outside 
the light cylinder, and this has to do with the global topol- 
ogy of the reconnecting magnetic field. Equations [2] imply 
that, next to the region of reconnection in the striped 
wind at large distances (r :$> ri,), it is the toroidal mag- 
netic field component that is the dominant one. Several 
authors, therefore, opted to ignore the poloidal magnetic 
field component altogether and instead study reconnec- 
tion as if it were taking place exclusively in the toroidal 
magnetic field component. It is as if the star 'shoots out' 
rings of a toroidal magnetic field with different 0-direction 
depending on whether they originate in the north or south 



loannis Contopoulos: Magnetic topology 



3 



magnetic hemispheres. The problem with that approach 
is that, by ignoring the global topology of the poloidal 
magnetic field component, we miss the crucial fact that 
the source of the toroidal magnetic field component is the 
azimuthal winding of the poloidal one. 

This is not without significance. All field lines orig- 
inate on the surface of the star, and there is a finite 
amount of magnetic flux that crosses the light cylinder and 
would be stretched out to infinity, were it not for the pres- 
ence of equatorial current sheet reconnection. Obviously, 
as toroidal magnetic field annihilation proceeds, less and 
less poloidal magnetic flux extends to larger distance^ 
As a result, less and less open field lines are available for 
reconnection or for carrying with them the stellar spin- 
down energy in the form of Poynting flux. We may thus 
generalize Eqs. [T] & [2] and introduce a simple power law 
scaling 



(3) 



«'(r,6') « -cos6l) — 



where e ^ 1 is a small positive parameter that obviously 
depends on the (yet unknown) physics of gradual recon- 
nection in the equatorial sheet. At this stage, we can only 
take it to be a free parameter. Equation[3]is an ad hoc gen- 
eralization of Eq. [1] and is not based on current theory or 
observations. It may, however, be valid as a toy model that 
describes toroidal magnetic fleld annihilation and poloidal 
magnetic field reconnection in the vicinity of the light 
cylinder. Obviously, in order to study how fast toroidal 
magnetic field annihilation proceeds with distance, one 
needs to solve a formidable non-force-free (beyond the 
fast magnetosonic surface) relativistic MHD problem that 
incorporates the microphysics of field annihilation. It is 
interesting to note that recent time-dependent numerical 
simulations (jBucciantini et al. 20061) suggest that recon- 
nection may be sporadic, spontaneous, and limited in the 
vicinity of the closed line zone. 

Having acknowledged the limitations of our toy model, 
it is now straightforward to show that 



|S.(r,( 



Bl 



We conclude this section by emphasizing that an- 
nihilation of the toroidal magnetic field in the equa- 
torial current sheet may be studied only in conjunc- 
tion with studying the poloidal structure of the mag- 
netospheric field. To our knowledge, this has not been 
taken into account in any existing model of extended 
equatorial field annihilation (e.g. Lyubarsky fc Kirk 200 1[ 
[Romanova, Chulsky fc Lovelace 2005 ). 



3. Evolution of the pulsar magnetosphere with 
neutron star spindown 

As we argued in Contopoulos (2007), the rate of equato- 
rial current-sheet reconnection in the vicinity of the light 
cylinder may be indirectly measurable through observa- 
tions of the neutron star's braking index 



(5) 



In that work, we dealt with non-ideal MHD regions 
through the notion of the effective magnetic diffusivity. In 
the present work, we give a simpler picture without men- 
tioning the magnetic diffusivity in the equatorial current 
sheet. 

In the above picture of finite plasma conductivity, mag- 
netic field lines enter the equatorial current sheet at a 
grazing angle, i.e. \Br\ ^ B^ outside the current sheet. 
Electromagnetic energy is continuously 'pumped' into it 
and is lost in accelerating particles in the current sheet 
(see next section). In other words, equatorial toroidal mag- 
netic field annihilation prevents poloidal magnetic flux 
from reaching 'infinite' distances. The picture we describe 
here may even be thought to take place in steady state 
were it not for the fact that the neutron star actually 
spins down. 

Why is the neutron star spindown so important? We 
remind the reader that, inside the light cylinder, the pulsar 
magnetosphere consists of a corrotating region of closed 
field lines and of a certain amount 5'opcn of open poloidal 
magnetic flux which extends beyond the light cylinder. It 
_^^s well known that the spindown rate tl depends on f2 and 



open 



as 



\B^{r,e)\ - BLsin( 



(4) 



(6) 



Note that Bg is everywhere of the same sign. This com- 
ponent of the magnetic field, which up to now has been 
considered insignificant, may have interesting observable 
implications. In fact. lPetri fc Kirk 20051 present a plausi- 
ble physical model for the Crab pulsar, where high energy 
radiation originates in the equatorial current sheet at a 
distance of about lOr^. In their model, they require a 
magnetic field component Eg roughly 10% of S,/,, which 
is not very different from what Eqs. 2] yield for e<l. 



In other words, toroidal magnetic field annihilation pro- 
ceeds in the north-south direction. 



( [Contopoulos 2005| . Let us backtrack for a moment and 
assume that reconnection does not take place anywhere 
in the pulsar magnetosphere. In that the neu- 

tron star spins down, the light cylinder moves out at a 
certain rate (roughly 1 km per year), but the region of 
closed field lines cannot grow, since the absence of recon- 
nection prevents open field lines from closing and accu- 
mulating on top of the closed line region, and thus Eq. [T] 
implies that 5'opon = = const. According to Eqs. [5]fc 
[51 the neutron star braking index would be measured to be 
equal to 1, which is not what is observed. Now let us once 
again consider the physically more interesting case where 
reconnection does take place in the vicinity of the light 
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cylinder, only it is confined inside the equatorial current 
sheet. In that case, as the neutron star spins down and the 
light cylinder moves out, it encompasses more and more 
poloidal magnetic flux because, as we argued, in the vicin- 
ity of the light cylinder magnetic flux closes with distance 
at a rate given by Eq. [3l When the light cylinder moves 
out from TL io rL + 5rL, a certain amount of poloidal 
magnetic flux. 



(5*opc 



5VL 

77 



(7) 



is now found contained inside the new light cylinder. 
Under ideal MHD conditions, however, the stellar spin- 
down energy is transferred only along field lines that cross 
the light cylindeJl; therefore, our result in Eq. [7] modifies 
the neutron star spindown rate by an amount 



5h = {l + 2e)^6n , 

which yields a braking index value of 
n = 1 + 2e . 



(8) 
(9) 



The reader may have noticed that our present discus- 
sion differs from the one in Contopoulos (2007), where 
we introduced an effective magnetic diffusivity 77, which 
is responsible for both the equatorial reconnection of the 
toroidal magnetic field and the inward diffusion of the 
poloidal magnetic field lines. In the present work, we again 
consider magnetic field energy losses in the equatorial cur- 
rent sheet, but we also argue that the growth of the closed- 
line region is directly related to the advancement of the 
light cylinder to larger distances as the neutron star spins 
down. 

We repeat once again that toroidal magnetic field anni- 
hilation is related to how fast poloidal magnetic field lines 
reconnect as we move away from the light cylinder (Eqs.[3] 
& IH). This is very important, because the efhciency of 
poloidal magnetic field closure is something that may be 
measurable through the pulsar braking index, and thus 
yield an estimate of the physical processes that lead to 
magnetic field reconnection in the pulsar magnetosphere. 

4. Particle acceleration in the equatorial current 
sheet 

Equations H] are valid outside the current sheet disconti- 
nuity. Inside the sheet, Br and change sign, meaning a 
neutral layer develops where the only remaining magnetic 
field component is Bg. In other words, when we generalize 
Bogovalov (1999), reconnection may take place at various 
positions along the current sheet discontinuity. All field 
lines originate on the surface of the star; therefore, as we 



^ As we show in the next section, a radial electric field equal 
to Er = Be{r /tl) appears inside the equatorial current sheet. 
Obviously, inside the light cylinder, E _L B and Er < B, so 
electromagnetic field energy cannot be dissipated through di- 
rect particle acceleration in the equatorial current sheet. 



argued, they carry with them the information that the 
neutron star is rotating at angular velocity 17 throughout 
the region where ideal MHD is a valid physical approxi- 
mation. In axisymmetry, a poloidal electric field develops 
perpendicular to the poloidal magnetic field because of the 
neutron star rotation, with magnitude 



E 



rfl sin 9 



[Bt+B, 



(10) 



fe.g. IBeskin 1997]) . Magnetic flux surfaces are, therefore, 
isopotential surfaces, with potential difference between 
them given by 



5V 



2tt'c 



(5* 



(11) 



where 5'^ is the magnetic flux contained between the two 
fltix surfaces. This is valid along each field line, all the 
way from its footpoint on the stellar surface to the point 
where that field line reaches the equatorial current sheet 
and reconnects (if that particular one does) with the cor- 
responding field line from the opposite magnetic hemi- 
sphere. According to our discussion above, this implies 
that the magnetospheric potential drop manifests itself 
along the equatorial sheet wherever poloidal field lines re- 
connect. Another way to understand this is that, wherever 
reconnection takes place, a radial electric field equal to 



Er = '-^Bs 
rL 



(12) 



develops in the interior of the current sheet where B^ ^ 
Bff, ~ 0. This is the reconnection electric field compo- 
nent mentioned in Kirk (2004) and Romanova, Chulsky & 
Lovelace (2005) that accelerates particles in the current 
sheet. The new result in our present work is that Er is 
directly proportional to Bg ^ Q in the current sheet. 

As long as a particle stays in the interior of a recon- 
nection region in the equatorial current sheet, it will be 
accelerated to an energy equal to eSV , where e is the elec- 
tron charge, and 5V is the potential drop associated with 
the poloidal field lines at the two ends of the accelera- 
tion region. Current sheets are known to be unstable, so 
we expect that the pulsar equatorial current sheet will 
also be fragmented into several short acceleration regions. 
Nevertheless, one may argue that some particles may man- 
age to profit from several of these acceleration regions. In 
such a stochastic acceleration, some particles may man- 
age to get accelerated to the maximum energy available 
in the open pulsar magnetosphere, i.e. to the energy due 
to the potential drop across field lines that originate on 
the magnetic pole and the edge of the neutron star polar 
cap, namely ~ 10^^ eV \De Jager fc Harding 1992D . 

We argue that synchrotron radiation reaction does not 
inhibit the acceleration. We performed a straightforward 
calculation of relativistic electron/positron orbits as they 
enter the current sheet (Figs. 1 & 2) and realized that 
their orbit becomes more and more confined around the 
center of the current sheet where Br ~ B^j, ^ 0. Therefore, 
even if one is willing to keep our toy model scaling all the 



loannis Contopoulos: Magnetic topology 



5 



way from to r 3> r/,, the synchrotron radiation reaction 
due to the Bg component of the magnetic field will limit 
the acceleration only when the particles reach a Lorentz 
factor on the order of 



10^ 



1/2 



B 



-1/2 



> 10^ , 



where, B^ is measured in Gauss. Thus, synchrotron radi- 
ation is unimportant inside the pulsar wind termination 
shock. 

We end our discussion with an order-of-magnitude 
estimate of the acceleration length /acc in the vicinity 
of the light cylinder. Electromagnetic energy flux en- 
ters the current sheet from above and below at a rate 
of 2{ErBfj,/4:Trc){Bg/Br)- This is absorbed by the par- 
ticle flux 2KnQjc{Bg/ Br), which also enters the current 
sheet from above and below. Here, k ^ 10"^ — 10^ is the 
multiplicity of pair production in magnetospheric gaps 
( Daugherty fc Harding 1982 IGurevich & Istomi n 1985). 
and ncj ~ flBr/ (47rce) is the Goldreich-Julian number 
density ("Go ldreich &: Julian 1969| . If we assume an accel- 
eration length Zacc ^ a straightforward energy balance 
in the current sheet yields 



2Ke 



(14) 



A (typical) value of e on the order of one half yields ^acc ^ 
lQ~*r ^ r. On the other hand, if we know Zacc from a nu- 
merical simulation of the pulsar current sheet fragmenta- 
tion (e.g. lCrew et al 20011 Bessho fc Bhattacharjee 2005 ), 
Eq. [14] yields e, which determines the global topology of 
the pulsar magnetosphere for various values of the multi- 
plicity parameter k. We thus predict that the smaller the 
efficiency of pair production, the larger e is. Unfortunately, 
the 6 pulsars with the best-measured values of the braking 
index are all strong pair creators, and therefore, our pre- 
diction cannot be tested observationally using only those 
particular 6 pulsars. 



5. Summary 

In the present work, we argue that toroidal magnetic field 
annihilation in the equatorial current sheet of the pul- 
sar magnetosphere cannot be studied independently of the 
poloidal magnetic field. This association allows us to re- 
late the pulsar braking index to the rate at which poloidal 
magnetic field lines close with distance in the vicinity of 
the light cylinder. This also allows us to estimate the re- 
connection electric field that develops in the interior of the 
equatorial current sheet. We show that this electric field 
accelerates electrons and positron in opposite directions, 
and we predict that there will be little or no synchrotron 
radiation in the region upstream of the wind termination 
shock. In fact, some particles may even profit from the 
full potential drop in the open line region and thus reach 
energies on the order of 10^^ eV. Finally, we associate the 
efficiency of equatorial reconnection to the efficiency of 
pair production in the pulsar magnetosphere. 
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Fig. 1. Calculated electron orbits that enter the equatorial 
current sheet 2; = from above. The units are in centime- 
ters. Wc consider only the axisymmetric case ■& = 0. The 
injection region is at r = IOtl- Here, Be/ Br = 10-^ on the 
surface of the current sheet. We took tl = 1500km and 
Bl ~ lO^G (values representative of the Crab pulsar). 
Initially, the relativistic electrons move together with the 
magnetic field with Lorentz factor equal to ~ 100. The 
Lorentz factor grows linearly with the acceleration dis- 
tance. We observe that the particles do not escape from 
the current sheet (unless of course the current sheet is 
terminated), and are confined more and more around the 
center of the current sheet where B = Be9. 



Fig. 2. Same as Fig. 1 for positrons. 



